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ABSTRACT 

Context. Herbig Ae stars (HAe) are the precursors of Vega-type systems and, therefore, crucial objects in planet formation studies. 
Thus far, only a few disks associated with HAe stars have been studied using millimetre interferometers. 
Aims. Our aim is to determine the dust evolution and the lifetime of the disks associated with Herbig Ae stars. 

Methods. We imaged the continuum emission at ~3 mm and ~1.3 mm of the Herbig Ae/Be stars BD+61154, RR Tau, VY Mon and 
LltHc 198 using the Plateau de Bure Interferometer (PdBI). These stars are in the upper end of the stellar mass range of the Herbig 
Ae stars (M„>3 M Q ). Our measurements were used to complete the Spectral Energy Distribution (SED). The modelling of the SED, 
in particular the FIR-mm part, allow us to determine the masses and dust properties of these disks. 

Results. We detected the disks associated with BD+61 154, RR Tau and VY Mon with disk masses of 0.35 Ms, 0.05 M Q and 0.40 M Q 
respectively. The disk around LkHo' 198 was not detected with an upper limit to the disk mass of 0.004 M Q . We detected, however, 
the disks associated with the younger stellar objects LkHo- 198-IR and LkHo- 198-mm that are located in the vicinity of LkHQ- 198. 
The fitting of the mm part of the SED reveal that the grains in the mid-plane of the disks around BD+61 154, RR Tau and VY Mon 
have sizes of ~ 1-1000 pm. Therefore, grains have not grown to centimetre sizes in these disks yet. 

Conclusions. These massive (M„>3 M Q ) and young (~1 Myr) HAe stars are surrounded by massive (>0.04 M G ) disks with grains 
of micron-millimetre sizes. Although grain growth is proceeding in these disks, their evolutionary stage is prior to the formation of 
planetesimals. These disks are less evolved than those detected around T Tauri and Herbig Be stars. 
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1. Introduction 

HAeBe stars, pre-main sequence intermediate mass objects 
(M»~2-8 Mo), share many characteristics with high-mass stars, 
such as clustering, presence of photo dissociation regions (PDR), 
etc. HAeBe stars are, however, much closer to us and less em- 
bedded than higher mass stars, which makes the detection of cir- 
cumstellar disks around these stars much more feasible, and key 
for the understanding of the mechanisms of massive star for- 
mation in general. In addition, determining the lifetimes and 
disk dispersal mechanisms of circumstellar disks is important 
for planet formation studies. Disk lifetimes set an upper limit 
on the time available for the assembly of planetary systems. 
Determining the details of the gas and dust dispersal is neces- 
sary to establish the inital conditions for planet formation. 

We are carrying out a comprehensive search for circum- 
stellar disks around HAeBe stars using the IRAM Plateau de 
Bure Interferometer (PdBI). Our final goal is to investigate the 
properties of the disks around intermediate mass stars to even- 
tually determine their occurrence, lifetime, and evolution. We 
have completed the observations toward 1 1 HAeBe stars. The 
results co rresponding to the first 7 s t ars have already be e n pub - 
lished by Fuen te et all d2003l 120061) . lAlonso-Albi et al. (2008), 



and lAlonso-Albi et al.l J2009) which will be referred to hereafter 
as Paper I. Four stars are presented in the present paper. 

The first 7-star sample was mainly composed of Be stars in 
the upper end of the intermediate-mass stellar range (M»>7 M Q ). 
Our results showed that the disks associated with these hot stars 
have masses a factor of 5-10 lower than the masses of the disks 
associated with Tauri (TT) stars. Since all these massive stars are 
very young, <1 Myr, we concluded that disk dispersal is signif- 
icantly more rapid at higher stellar masses (>7 M ) than in the 
lower mass stars. This result is also consistent with recent con- 
clusions by iBerne et ail (120091) based on PAHs emission from 
disks. They found that the PAH emission observed toward hot 
HBe stars (>7 M Q ) is associated with the nebula instead of the 
disk proving either the non-existence or the existence of smaller 
disks around these stars. Regarding the dust properties, our re- 
sults proved the existence of large grain sizes (~1 cm) in the 
midplane of the disks associated with these massive, >7 M , 
stars. If planet formation occurs in these stars, it should occur in 
timescales of a few 10 5 yr. 

Several groups have been obser ving the disks associated with 
Herbig Ae stars in the last decade ([Testi et al.ll2003t iPietu et al.l 
20031 120051 l2006t lHamidouche et alJl2006t lHamidouchd| 201( ' 



Isell a et alj20()7t I Al onso-Al bi et alj2008tlChapillon et al]2 008; 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



Banzatti et alJ 120111) . These stars have typical ages of a few 
Myr and stellar masses between 2 and 3 M Q . Evidences for 
grain growth have been found in many of these disks (see 
e.g. iTesti et al.ll200l llsella et alJl2007t lAlonso-Albi e"t~aT1l2008l: 



2 



J. Boissier et al.: Massive young disks around Herbig Ae stars 



Table 2. Plateau de Bure observations 
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m Tv/hpam 
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bD+61 154 


1.1 


1.12x0.7 


Dec. 2008 


0.5 




3.0 


2.68x2.54 


Apr. 2008 


0.3 


LkHa 198 


1.3 


1.56x0.97 


Mar. 2009 


0.5 




3.5 


1.40x1.11 


Jan. 2009 


0.1 


RRTau 


1.3 


2.27x1.75 


Apr. 2008 


0.8 




3.5 


1.35x0.74 


Feb. 2009 


0.1 


VYMon 


1.4 


1.98x1.89 


Oct. 2009 


0.6 




3.4 


7.50x4.45 


Aug. 2009 


0.4 



_ Contours: 2ct 



Banz atti et al.l l20lTh . Thus far, there is no evidence for a dif- 
ferentiated evolution from that of the disks around T Tauri stars 
(TTs). No correlation between the disk mass and neither the stel- 
lar mass nor the stellar age have been found in TTs a nd Herbig 
Ae stars (see e.g. lNatta et al.ll2007t iRicci et al . 2010bE). 

These observat ional results agree with recent photo- 
evaporation models iGorti & Hollenbachl J2009) which predict 
that the disk lifetimes are similar for disks associated to stars 
with stellar masses in the range of 0.3-7 M , but rapidly decline 
for stars of higher masses (Be stars). However, observational 
data are still very scarce, specially in the -3-7 M Q stellar mass 
range where the change between one regime and the other oc- 
curs. We present continuum interferometric observations using 
the PdBI and the subsequent SED modelling of four sources in 
this stellar mass range. Our goal is to get some insights into the 
physical characteristics of the circumstellar disks around these 
stars that are in the borderline between Herbig Ae and Be stars. 

2. Observations 

We carried out interferometric observations in the continuum at 
~ 1 .3 mm and -3.0 mm toward the 4 stars listed in TableQ] These 
objects are HAeBe stars with stellar masses of ~3-7 M located 
at a distance of 600-800 pc from the Sun. All the sources were 
observed in 2008-2009 using the PdBI in different configura- 
tions. A summary of the observations is presented in Table |2] 
including dates, wavelengths and beam sizes. The data reduc- 
tion and analysis were performed using the GILDAS software 
package. The continuum maps presented in the following were 
built summing all the channels free of line emission in the data. 
Natural weighting has been applied to the measured visibilities. 

3. Results 

The observed ~1.3 mm and ~3.0 mm continuum images are 
shown in Fig. Q] We detected 3 out of the 4 targets observed 
with the PdBI. We measured the fluxes by fitting different source 
models to the visibilities in the Fourier plane (wv-plane). We se- 
lected the source models according to the signal to noise ratio 
(S/N) of the observations : point source for lower S/N, Gaussian 
distribution (circular or elliptical depending on the source) for 
higher S/N values. All the fit results (positions, fluxes and width 
for Gaussian fits) are reported in Table [3] 

BD+61 154 was previously observed with OVRO at 2.7 mm 
by Mannings & Sargent who measured a flux of 1 1 .2 mJy (t hese 
observations were quoted in the review bv lNatta et al.S2 000). We 
measured a flux of 1.7+0.3 mJy at 3.0 mm using the PdBI. This 
large difference in the flux is probably due to the different fre- 
quencies (a factor of 1.4 assuming a spectral index of 3) and to 
the higher angular resolution of our observations that filter out 
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Fig. 1. Interferometric maps of the continuum emission of 
BD+61154, LkH al98, RR Tau and VY Mon around ~3 mm 
(left panels) and ~1.3 mm (right panels) as observed with 
the Plateau de Bure interferometer. The exact observing wave- 
lengths are indicated in the lower right corners. The level spac- 
ings (in units of map noise) are indicated in the upper left corners 
of the frames and were chosen according the signal to noise ra- 
tios of the detections (see Table|2]for the noise values). The syn- 
thesized beams are presented in the bottom left corners of the 
panels. The beam sizes and observing dates are given in Table|2] 
The measured fluxes (results of fits in the Fourier plane) are re- 
ported in Table [3] 



the extended emission. IWendkerl (l 995 ) published a 1.3 mm con- 
tinuum flux of 60 mJy based on single-dish observations carried 
out with the 30m telescope in this same source. Our interfero- 
metric observations recover ~50% of the continuum single-dish 
flux. These young stars are still embedded in the parent core and 
the single-dish emission is frequently dominated by the emis- 
sion from the envelope (see Paper I). According to our fit of the 
~1.3 mm data, the size of the emitting region is ~260 AU. 
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Table 1. List of sources 



Star 


Sp. 


RA(J2000) 


Dec(J2000) 


d 


Star Mass 


Age 1 






h:m:s 




pc 


Me 


Myr 


BD+61154 


B8 


00:43:18.257 


+61:54:40.13" 


650 


4.0 


<1 2 


LkHff 198 


B9/A5 


00:11:25.970 


+58:49:29.10" 


600 


3.5 


1.5 


RRTau 


A2 


05:39:30.502 


+26:22:26.90" 


800 


3.6 


1.3 


VY Mon 


B8 


06:31:06.937 


+ 10:26:05.04" 


800 


6.9 


0.2 



1 Stellar ages calculated using the Siess tracks and the stellar data in Table A. 13 of lAlonscT Albietal.l f2009). 

2 See discussion in Sect 5.1. 



We did not detect continuum emission toward LkHa 198 
with an rms of 0.1 mJy/beam at 3.5 mm and 0.5 mJy/beam 
at 1.3 mm. Our upper limit to the -3 mm continuum flux is 
a facto r of 20 lower than that obtained by Idi Francesco et al.l 
dl997l) . Within the field of view of the LkHor 198 image, 
we detected two millimetre peaks that correspon d to the mid- 
IR sou rce LkHa 198-IR previously detected by Lagag e"et al.l 
d 1993b and the submillimetre protostar LkHa 198-mm detected 
bv lSandell & Weintraubl d!994l) . LkHQ-198-mm is located at the 
border of the primary beam at 1.3 mm and for this reason, the 
flux is quite uncertain. 

Our interferometric measurements toward RR Tau are the 
only ones at millimetre wavelengths. We detected a flux of 
17.9 mJy at 1 .3 mm and flux of 2. 1 mJy at 3.5 mm and estimated 
the size of the emission region to -720 AU from the 1.3 mm 
data. 



iHenning et al.l d!998f) measured a flux of 112 mJy for the 
dense core in VY Mon using the IRAM 30 m telescope. Our 
1 .3 mm PdBI observation gives a flux of 60.4+ 1 .4 mJy, therefore 
our interferometric observations filtered out about 50% of the 
core flux. Our 1.3 mm data suggest a size of -980 AU for the 
emitting region in this source. 

All the sizes derived above (-260 AU for BD+61154, 
-720 AU for RR Tau, and -980 AU for VY Mon) are typical 
for circumstellar disks around intermediate mass stars. Although 
some emission from the envelope can be detected with our inter- 
ferometric observations, the small source sizes confirm that the 
detected emission arises mainly from the disk. This interpreta- 
tion is also reinforced by the fact that the SEDs we have con- 
structed with our interferometric fluxes are well fitted assuming 
that all the emission detected in interferometry arises from the 
disk (envelope emission is negligible). 

A fraction of the interferometric flux measured at -1.3 mm 
and -3.0 mm could be due to the stellar wind emission in- 
stead of the dust thermal emission from the disk. We have es- 
timated the contribution of the free-free emission at millimetre 
wavelengths by extrapolating previous observations at cen time- 
ter w avelengths assuming the power-law F v <x y +0 - 6 (e.g. lOlnonl 
1 1975b . The spectral index, a - +0.6, is expected for an opti- 
cally thick isotropic wind. This value is fully consistent with 
those measured in Herbig Be s tars using VLA observations by 
lAlonso-Albi et al.l (l2008l l2009t) . I t is also consis t ent, within the 
errors, with the values derived by ISkinner et al .1 (119931) in most 
HAeBe stars, specially towards the most massive stars of their 
sample. Only in two cases, V645 Cyg and TY Cra, the spec- 
tral index clearly differs from this value. The emission towards 
V645 Cyg was over-resolved and the derived spectal index was 
a — 0.2+0.3. This value is consistent with optically thin free- 
free emission or the emission from an HII region. In the case 
of TY Cra, a negative spectral index of -1.2 is found as well 
as signs of vari ability pointing to a non-thermal emission mech- 
anism (see also iLommen et al.1 120091) . iNatta etafl d2004l) found 



a - +0.6 in HD 163296 and a value of -0.0 for HD 35187. We 
consider that a - +0.6 is a reasonable assumption to estimate 
the possible contribution of the gas continuum emission at mil- 
limetre wavelengths. 

In Table [3j we show the millimetre spectral index derived 
from our data after subtracting the free-free emission, when its 
contribution is not negligible. The millimetre spectral index is a 
good indicator of the grain sizes in the disk midplane. The val- 
ues in Table [3] are larger than those found in the massive stars 
reported in Paper I, suggesting that the grain sizes are lower: we 
are dealing with a class of less evolved disks. 



4. Modelling the SED 

We constructed the whole SED toward our sample by combin- 
ing our PdB data with other results available in the literature 
and Spitzer photometry fluxes taken in the public archive. These 
source fluxes were extracted in a circular aperture of 10" radius. 
Since these stars are bright, we took an off-source flux far away 
(~ 1') from the source to avoid contamination. The off, contain- 
ing zodiacal light and ISM contamination is subtracted to the 
on-source flux. The accuracy is about 10% for this method. 

At optical and NIR wavelengths, observations are heavily af- 
fected by extinction. To correct this reddening we subtract the 
standard B-V col or correspondin g to the spectral type of each 
star, as given bv iJohnsonl (fl966). to the observed B-V color. 
This color excess E is used to calculate the visual extinction 



A,, = R v E assuming a rati o of total to selective absorption R v 
= 3.1 (Cardell i et al.lll989l) . With these two parameters we es- 
timate the extinction corresponding to a particular observatio n 



timate the extinction corresponding to a particular observation 
using the parameterized extinction law of Cardelli et al. (1989), 
which can be applied for wavelengths between 0.12 and 3.5 
fim. Most HAeBe stars are well known variable stars, a prob- 
lem that becomes relevant in our case since we are using pho- 
tometry recovered from the literature in the last decades. The 
variability affects the B-V color and produces some additional 
uncertainty in the extinction estimates. The disks have been fit- 
ted using the passive irradia ted circumstellar disk model de- 
scribed by lAlonso-Albi et al.l (120081) that allows us to consider 
different grain populations for the midplane and surface layers. 
Grain populations are characterized by the silicate/graphite mix- 
ture, the maximum grain size (a max ), and the slope of the grain 
size distribution («). A standard grain mixture (SM: 86% silicate 
and 14% graphite) and a value of n — 3.5 are assumed for the 
midplane and surface layers. We allowed a max to vary between 
0.1 fim and 1 cm in the midplane layer. Observations at mil- 
limetre and centimetre wavelengths are not sensitive to grains 
of sizes larger than a few centimetres. During the sedimentation 
process, large grains (>1 mm) are expected to accumulate in the 
midplane while the surface layer remains populated by smaller 
grains. For this reason, in the surface layer we only consider the 
values a max =0.1, 1, 10, 100 fim. The gas to dust ratio is assumed 
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Table 3. Results of the fits to the visibilities. Uncertainties on the output parameters are given in parenthesis. 





Thi itipti on 


\ 


Source Position 


HPW 


PA 


Flux 








mm 


RA(J2000) 


Dec(J2000) 






mJy 




BD+61154 


E-GAUSS 


1.1 


00:43:18.257(0.001) 


61:54:40.13(0.01) 


0.39(0.06)x0. 19(0.06) 


27(13) 


33(1) 


2.9(0.6) 




POINT 


3.0 


00:43:18.160(0.028) 


61:54:39.90(0.19) 






1.7(0.3) 




RRTau 


C-GAUSS 


1.3 


05:39:30.516(0.004) 


26:22:26.83(0.05) 


0.9(0.2) 




17.9(1.7) 


2.2(0.4) 




C-GAUSS 


3.5 


05:39:30.514(0.003) 


26:22:26.94(0.04) 


0.3(0.1) 




2.1(0.2) 




VYMon 


E-GAUSS 


1.4 


06:31:06.938(0.001) 


10:26:04.61(0.01) 


1.22(0.06)x0.77(0.07) 


-21(6) 


60.4(1.4) 


4.0(0.2) 2 




POINT 


3.4 


06:31:06.926(0.003) 


10:26:04.72(0.04) 






9.7(0.2) 





' The fit of the Fourier Transform of point source in the t(v-plane has 3 free parameters: the source RA and Dec positions relative to the phase 
center and its flux. Fitting a Gaussian distribution has additional parameters: the width (HPW) for circular Gaussian (C-GAUSS) and the major 
and minor axis for an elliptical Gaussian (E-GAUSS), as well as its orientation with respect to the north (Position angle PA). 
2 This value has been derived after subtracting the contribution of the free-free emission (see text, Sect. 15.3) . 



to be 100. In addition to a max , we vary the inclination angle of 
the disk (i, in the range 10-80°), the rim temperature (Tnm, in the 
range 1200-2000 K, the sublimation temperatures of silicates 
and graphite being 1500 and 2000 K, respectively) and the index 
of the density law (p, in the range -2.5-0). Finally, the disk outer 
radius (R out ) and dust mass (M^k) vary in ranges defined accord- 
ing to a first guess fit of the SED of the sources. These ranges 
are, respectively : 30-200 AU and 1.0xl0- 4 -10.0xl0- 4 M Q for 
BD+6154; 30-150 AU and 2x10 4 -32.0xl0 4 Mq for RRTau; 
100-500 AU and 2.0xl0- 4 -50.0xl0~ 4 M for VY Mon. 

The goodness of the fit is estimated with a dispersion param- 
eter (x) defined as: 



_ / 1 ^ ^ observed ^ modeled)^ 
V n r observed 



(1) 



where F is the flux and where the sum is extended for all de- 
tections in the SED from 3 /im and beyond, being n the number 
of them (upper limits are not considered). Assuming that flux 
uncertainties are dominated by calibration errors, which we as- 
sume to be constant all over the SED (and of the order of 10%), 
this dispersion is proportional to the standard x- The same dis- 
persion definition was used in Paper I. We report in Table |4]the 
best fit parameter set, which corresponds to the minimum value 
of for the three sources we detected with the PdBI. Table |4] 
contains as well the limits (Min and Max values) of the intervals 
within which x is lower than ~1.3 Xmm- They were calculated 
by varying each parameter separately while keeping all the other 
to their best fit values. 

As Mdisk and a max affect only the fluxes in the FIR and mm 
ranges, where we have only few measurements (sometimes only 
our two PdBI fluxes at ~1.3 mm and ~3.0 mm), their changes 
impact faintly the normalized dispersion as defined in equa- 
tionfTJ To enhance their impact, we estimated their Min and Max 
values considering only the restricted SED range of A > 100/im. 
Our fits are shown in Figs. |2]to |4] In the following we comment 
the individual sources and the uncertainties in our estimates. 

5. Individual sources 

5.1. BD+61154 

BD+61154 (also design ated as MWC 419 or V594 Cas) 
has spectral type B8 (iHerbigl 1 19601) . luminosity 330 L 



(Hillenbrand et al 



(Hillenbrand et al 



19921), and is at a d istance of 650 pc 



1992h! lTestietal.l(ll998h estimated an age of 
0.1 Myr assuming T P ff =11220 K and a luminosity of 300 L . 
Later. lAlonso- AM et al.ld2009l) derived an age of 1 .4 Myr adopt- 
ing a lower value of the luminosity of 220 L . This low value of 



the luminosity cannot fit the NIR part of the SED. We need a 
luminosity of ~300 L to fit the SED, which implies that the star 
is younger than 1 Myr (see Table [TJ. 

Based on prev ious VLA observations published by 
ISkinner et al ] (119931) . we estimated that the maximum contri- 
butions of the free-free emissions are 0.44 and 0.8 mJy at 3.0 
and 1 . 1 mm respectively, which is negligible in comparison with 
the fluxes measured at the PdBI. Consequently, the free-free 
emission is not taken into account in the fit of the SED for 
BD+61 154. In Fig. |2^ we show our best fit to the SED of this 
star. The bump appearing in the SED between 1 and 4 microns 
is well reproduced assuming a rim temperature of 1800 K. The 
agreement between our model and observations is quite good 
and allows us to constrain some disk parameters. The best fit is 
obtained with an inclination of 70°, an outer radius of 50 AU and 
a density distribution slope of -0.7. Our millimetric data con- 
strain the dust disk mass to 3.5 x 10~ 4 M and the dust maximum 
grain size to 100 fim in the midplane layer. Within the errors, the 
radius of the disk obtained by fitting the SED of BD+61154 is 
consistent with that derived from the fitting of the PdBI ~ 1 .3 mm 
continuum emission (~125 AU), which reinforces our interpre- 
tation that the interferometric millimetre fluxes originate in the 
disk. 

One has to keep in mind however that the outer radius is 
determined by the FIR part of the SED and we have very few flux 
measurements in this region (see Fig. |2j. Photometry in the 100- 
1000 //m wavelength range is necessary to have a more accurate 
estimate of the outer radius. 



5.2. RR Tau 

RR Ta u is a high l y variable HAeBe star also class i fied as UVOR- 
tvpe (lTheetal.1 ll994t lOudmaiier etafl 120011: iRodgers et all 
2002), making difficult the estimate of its stellar parameters. It 
has been classified between B8-A4 b y different author s, with 
the most recent estimate as A2 (Hernand ez et al.l [20041) . Even 
more unce rtain is its distance. ITesti et al.l (1 19981) gave a value 
of 800 pc. iBlondel & DUel d2006l) modeled th e UV spectra and 
derive d a distance of 600 pc. More recently, iMontes inos et al. 
(120091) modeled the SED at UV and NIR wavelengths and de- 
rived a distance of ~2100 pc. The authors discussed, however, 
that the obtained value is very dependent on the assumed ex- 
tinction and the distance could be lower if a circumstellar disk 
blocks the stellar light. In this paper, we have adopted a distance 
of 800 pc and T e yy=9000 K that seems a good compromise tak- 
ing into account the data in the literature. The values of disk 
mass and dust spectral index derived from our modelling do not 
strongly depend on on the stellar effective temperature providing 
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Fig. 2. a) The model fitted to the SED of BD+61154 is shown with all the model contributions: interior layer (blue medium sized 
dashes), surface layer (red short and medium sized dashes), inner rim (pink dotted line), star (green short and large sized dashes), 
and total sum of previous contributions (black solid line), b) This panel illustrates how the SED changes when one varies the value 
of the maximum grain size in the interior layer (a„ ra v ). The inset show the projected dispersion curve (see Eq. Q~]for the definition of 
the dispersion x) when only a max varies, all other parameters being fixed at their best fit values, c) The same as b), but varying the 
outer radius, d) The same as b), but varying the disk mass. The x axis of panels b) and d) has been reduced to A > 100 //m, were the 
SED is more affected by a change of the parameters a max and disk mass, respectively. 



it is between 8000 K and 10000 K as expected for a B6-A4 star. 
The distance, however, has an important impact on the disk mass 
estimate, as it scales as D 2 . 

Using a distance of 800 pc, we got a satisfactory fit of the 
SEP. Given the upp er limit on the continuum flux at 3.6 cm 
(ISkinner et al.f l 993) and the spectral index of +0.6 for the free- 
free emission, it could account for at most 0.3 mJy at 3.5 mm 
and 0.6 mJy at 1.3 mm, which would not change significantly 
the SED fitting process. Therefore, we did not include it in the 



fit. The bump appearing in the SED between 1 and 2 microns is 
well reproduced assuming a rim temperature of 1200 K. The ge- 
ometry of the disk is constrained by the mid infrared part of the 
SED and the best agreement between our model and the obser- 
vations is obtained for a disk outer radius of 80 AU, an inclina- 
tion of 70° and a dust density slope p=-1.5. Our interferometric 
measurements at millimetre wavelengths are the only ones that 
are not upper limits above 500 fim. They allow us to estimate 
the dust mass of the disk, Md ust =5 x 10~ 4 M and the maximum 




size of dust grains in the midplane, a max =1000 fjm. The mea- 
surements performed by IRAS and ISO around 200 fim are not 
used in our fit because given their large beams, the fluxes mea- 
sured in these observations are likely dominated by the cloud 
emission. Assuming a standard gas/dust ratio of 100, the total 
dust+gas mass of the disk is 0.05 M Q . The size of the disk as de- 
termined by fitting the SED is smaller than that derived from the 
millimetre interferometric observations. The large dispersion of 
the flux measurements between 100 fim and 1000 yum cast some 
doubts about the accuracy of this estimate. On the other hand, 
the disk could have a complex geometry in the outer part that 
it is not well described by our simple model. The existence of 
spiral arms and r ings are no t unusual in HAeBe stars (AB Aur: 
iPietu et aUl200l MWC297: lAlonso-Albi et aD2009h 



5.3. VYMon 

VY Mon is a B8 star located at a distance of 800 pc. It is the most 
massive star in our sample. We have not been able to find a good 
fit for A < 1 jjm probably because the poorly known extinction 
and photospheric fluxes in this variable star. This, however, does 
not affect our fit to the disk since the stellar effective temperature 
and the distance are well known. Th e free-free emission that was 
detected at 6 cm bv lWendkerl (1 19951) corresponds to a flux of 8.4 
mJy at 3 mm after extrapolating with a 0.6 spectral index. This 
flux dominates the emission we observed with the Plateau de 
Bure. At 1 mm, the free-free emission contributes t o 14.3 mJy, 
i.e. a third of the observed flux. lHenning et all (I 19981) imaged the 
region at 1.3 mm using the IRAM 30 m telescope. They derived 
the existence of an extended component of ~53"x49" and a total 
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Wavelength (/um) Wavelength (fim) 

Fig. 4. Same figure as Fig.[2]for VY Mon. Note that in this case we also consider the free-free emission component (grey large sized 
dashes). 



flux of 0.25 Jy, and a compact component, <13" with a flux of 
0.11 Jy. The ISO 200 fan flux and the SCUBA 850 /urn fluxes 
are also expected to be dominated by the emission from the ex- 
tended component. At wavelengths larger than 100 fim, we only 
consider the interferometric observations for our fit. Our mil- 
limetric data constrain the dust mass (0.004 M ) and the dust 
maximum grain size to 1 /im in the midplane and 0. 1 yum in the 
surface of the disk. Assuming a standard gas/dust ratio of 100, 
the total dust+gas mass of the disk is 0.40 M . The disk size de- 
rived from the SED fitting is consistent with that estimated from 
the 1.3 mm PdBI images. 



5.4. LkHa 198 

LkHa 198 (V* V 633 Cas, IRAS 000 8 7+5833) is a HAeBe star 
located at 600 pc (IChavarria-K.lll985l) . lHillenbrand et al.1 (Il992h 
classified LkHa 198 a s a A5 pre main-sequence star, while 
iHernandez et alj (120041) found a spectral type of B9+2.5. 

A strong CO bipolar outflow is associated with this re- 
gion, but it is no t centered on LkHa 198, b ut to the north- 
west (ICanto et al.1 Il984t iBallv & Ladal Il983h . A small clus- 
ter of intermediate-mass young stellar objects is surrounding 
LkHa 198. V376 Cas is another HAeBe star located about 40" 
north of LkHa 198. An infrared co mpanion (LkHa 198 -IR) was 
detected ~6" north of LkHa 198 bv lLagage et al.1 (Il993l) . A sub- 
millimetre sourc e (LkHa 198-mm) wa s ident if ied 19" northwest 
of LkHa 198 bv lSandell & Weintraubl (fl994l) . lHaiiar & Bastienl 
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Table 4. Disk Models 
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1 Maximum deviation allowed from the minimum value of x to estimate 
the Min and max values. 

2 We have not considered values lower than 30 AU for the outer radius 
of the disk around BD+61 154.. 

3 We have not considered grains with a„ MV < 0.1 fan. 



(2000) classified this source as an extreme Class I object. 
Extreme class I objects have been proposed by Ladaj (Il99ll) to 
refer to early class I objects that are seen through an edge-on 
disk, thus explaining why they cannot be observed at N IR wave- 
lengths. BIMA observations by Matthe ws et alJ d2007h revea led 



that the single-dish outflow detected by Canto et al. d 19841) is 



indeed the combination of the several bipolar outflows associ- 
ated with these stars. Evidences for the existence of bipolar out- 
flows associated with V376 Cas and LkHa 198 were found. One, 
even two outflows seemed to be associated with LkHa 198, the 
most prominent in the CO emission with its axis in the northeast- 
southwest direction. A strong north-south outflow is associated 
with LkHa 198-mm. Although CO emission is detected toward 
LkHa 198-IR, there is no evidence of outflow in this IR source. 

We did not detect millimetre emission toward LkHa 198 at 
~1.3 mm and ~3.0 mm Intense emission is found, however, as- 
sociated with LkHa 198-IR and LkHa 198-mm (see Fig.Q]and 
Table[5]l. We used our 3<x upper limit to the ~1.3 mm emission 
to derive an upper limit to the mass of the disk associated with 
this HAeBe star. This limit has been calculated assuming op- 
tical thin emission at 1.3 mm, that is the usual case in disks 
around TT and Herbig Ae/Be stars, an averag e disk temperature 
of ~40 K in agreement with the predictions of Na tta et all (2000) 
for a B7 star, a dust opacity k\s mm = 1.0 cm 2 g~' and a gas/dust 
ratio of 100. We obtained that the gas+dust mass of the disk as- 
sociated with LkHa 198 is lower than 0.004 M Q . Alternatively 
we can have a small optically thick disk. In this case, our upper 
limit gives an upper limit to the disk size but does not inform 
about its mass. Assuming that the disk average temperature is 
~40 K, the radius of the disk should be <9 AU to be consistent 
with our 3-cr upper limit to the 1.3 mm emission. LkHa 198-IR 
and LkHa 198-mm are two intermediate-mass young objects. 
LkHa 198-mm dominates the flux at 1 mm of the region and 
drives an energetic bipolar outflow. Therefore, the existence of a 



circumstellar disk is expected. LkHa 198-IR is supposed to be 
a more evolved object. The two sources have millimetre spectral 
indeces ~2.4, consistent with the emission arising from a disk 
and with the values measured for BD+61 154 and RR Tau. In 
order to derive the disk masses, we have assumed temperatures 
ranging from 15^-0 K for the dust, that are reasonable values for 
disk associated with HAe stars. We obtained that disk masses of 
0.07-0.26 M , 0.42-1 .42 M Q for LkHa 198-IR and LkHa 198- 
mm respectively. 



5.5. Uncertainties in the disk mass, maximum grain size and 
disk outer radius 

The assumptions made in the SED fitting process introduce some 
uncertainties in the values derived for the disk properties. 

The grain size distribution is usually taken to be n(a) = «o X 
a~ n for a < a max , where no is a normalization factor and a max is 
the maximum grain size. Experimental studies of fragmentation 
find that the value of n could vary betwee n 1 .9 for low- velocit y 
collisions, and 4 for catastrophic impacts (Davis & Rvanlfl9 90). 

iTanaka etail (1 19961) argued that the standard n — 3.5 power 
law is a very general result that depends only on the assumption 
that the fragmentation process is self-s imilar and that the c ol- 
lision rate varies as a 2 . Simulations by ITanaka et alJ d2005l) of 
particle growth in stratified protoplanetary disks give vertically 
integrated size distributions that at late times are well described 
by a power law n = 3.0 for a < a max , plus a population of much 
larger bodies. The 1.3 mm/2.7 mm spectral index is considered a 
good indicator of the grain properties. In the Ray leigh- Jeans re- 
gion of the spectrum, the dust emission is proportional to y (2+ #> 
where j3 is the opacity index that can be accurately derived by fit- 
ting the observed SED. The value of ft in the submillimetre and 
millimetr e ranges is an excellent indicator of the grain size dis- 
tribution (Draine 2006). There is some degeneracy between the 
index of the grain size distribution, n, and the value of the max- 
imum grain size, a max . An index of n — 2.5 means that a larger 
number of grains are closer to the maximum grain size than for a 
more typical value of 3.5. Assuming n = 2.5 we would reach the 
same values of ft with grains slightly smaller . However, as pre- 
viously discussed by sev eral authors (see e.g. lAlonso"- Albi et al. 
2008; iRicci et aill2010bl) . for reasonable values of n, grain with 
sizes larger than few millimetre are needed to have values of /3 
lower than 1, as this is the case for RR Tau. 

The dust composition in the midplane cannot be inferred 
from the SED. For this reason, we assumed the standard mix- 
ture. The silicate feature at ~9.8 prni provides, however, some 
information about the dust composition in the surface layer. In 
all of our sources, the silicate feature is present and can be fitted 
with the standard mixture. The absence of this feature, however, 
does not imply the absence of silicate grains but the lack of sili- 
cate grains at a temperature of ~800 K. The silicate grains could 
be either too far from the central source and thus not heated to 
temperatures sufficiently high to allow emission in the mid-IR, 
or these grains could be too large to be heated efficiently. The 
dust composition and maximum grain size in the midplane de- 
termine the value of dust emissivity at 1.3 mm, Ki.3 mm , and con- 
sequently has an influence on the derived dust mass. In Fig. 5 
we show the value of Ki.3„ lm as a function of a max for different 
values of n and different grain compositions. For values of a max 
in the range 1 mm- 1 cm, all grain compositions and grain size 
distributions result in values of Ki3, nm that are within a factor of 
2-3 in agreement with the canonial value, 1 g~' cm 2 . For smaller 
grains, Ki.3 mm has a stronger dependance on the grain composi- 
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Table 5. Interferometric fluxes and masses for LkHff^S field sources 
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RA(J2000) 


Dec(J2000) 


S(3.5 mm) 
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1>15K T rf =40K 
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Fig. 5. Dust opacity at 1.3 mm as a function of the maximum 
grain size, a max , for different values of the grain size distribution 
slope, n {top) and for different dust compositions (bottom). Note 
that the variation of Ki.3 mm because of the different grain size 
distributions and grain composition is at most of a factor of ~10. 

tion and size distribution. The variation reaches at most a factor 
of 10. 

The outer radii calculated from the SED tend to be smaller 
than those derived from the interferometric data fitting proce- 
dure. In the case of BD+61154 and VY Mon, the difference is 
within the modelling uncertainties and the results of our mod- 
els are consistent with the observations. In the case of RR Tau, 
the outer radius calculated from the SED is definitely smaller. 
In the SED modelling we assume a surface density power-law 
truncated at the outer radius. Several authors have proposed a 
more realistic surface density distribution with a tapered ex- 
ponential edge. This kind of surface density distribution gives 
smaller outer radii and there fore less consistent with our obser- 
vations (Hughes et al. 2008). As discussed in Sect. 5.2, the large 
size measured by our interferometric observations in RR Tau is 
likely due to a more complex geometry of the outer part of the 
disk, with the possible contribution of a tiny envelope or some 
other extended structure such a secondary ring or spiral arms. 
Resolved images are required for a definitive answer. 

6. Discussion 

We present interferometric data on 4 HAeBe stars with stellar 
masses in the range 3-7 M Q and ages ~ 1 Myr. All of them except 
LkHa 198 have been detected and our SED modelling revealed 
that they are surrounded by massive disks (masses of about 0.03- 
0.4 Mq). Grains in these disks have grown only to sizes of ~ /um- 
mm. Although grain growth is proceeding in these disks, their 



evolutionary stage is prior to the formation of planetesimals. One 
possibility to explain the lack of disk detection in LkHa 198 is 
that t he age of th i s star is wrong because of the uncertain spectral 
type. Testi et al] d 1 998b estimated an age of 10 Myr assuming a 
spectral type of A4. If the star is older, most of the disk material 
could have been already dispersed by the star. Another possibil- 
ity is that the differences are due to the vicinity of young pro- 
tostars. The UV radiation and bipolar outflows associated with 
these YSOs could have truncated the disk. 

6.1. Correlation of the disk parameters with the stellar 
parameters. 

To investigate the possible dependence of the disk mass on the 
stellar mass and age, we updated our previous compilation of 
disks (Table A. 13 in Paper I). It is important to calculate the 
disk masses in a uniform way. Similarly to Paper I, we cal- 
culated all disk masses, including our own sample, by assum- 
ing optically thin emission at millimetre wavelengths, an aver- 
age disk tempera ture that depends on the stellar spectral type 
dNatta et al. 2000), a dust opacity k\ .3 mm =1.0 cm 2 g 1 , and a dust 
spectral index /J=1.0 in case we do not have ~1.3 mm obser- 
vations. We subtracted the free-free contrib ution in the sources 
previously studied by our team (thi s work; iFuente et al.l |2003 , 
l2006llAlonso-Albi et al.ll2008ll2009h 

In Fig.|6^, we represent the total disk masses (gas+dust) as a 
function of the stellar mass assuming a gas/dust ratio of 100. We 
only show in the plot the detections and the meaningful upper 
limits (disk masses <0.04 M ). For stars with masses of <3 M Q , 
there is a cloud of points around a mean value of ~0.04 M . Stars 
with masses of 3-7 M usually have massive disks (>0.04 M Q ). 
But there are three stars without disk detections, and with upper 
limits to the disk masses of <0.01 M . In this range of stellar 
masses, most of the detections come from this work and corre- 
spond to young stars. This could bias our statistics to larger disk 
masses. We need to complete our sample with older Ae stars to 
have a more realistic view of the disk evolution and determine 
the disk lifetimes in this stellar mass range. Most massive stars 
(>7 M ) have usually disk masses <0.04 M , consistent with a 
more rapid disk dispersal as we concluded in Paper I. 

It seems established that disk masses depend on the stellar 
mass because it determines the photoevaporation efficiency and 
consequently, the disk lifetime. Disk masses should also depend 
on the the stellar age because the star is expected to progres- 
sively disperse the gas and dust in the disk. In Fig. [6J5 we plot 
the total disk masses as a function of the stellar ages. Massive 
stars (>7 M ) are not included in this plot because it is not pos- 
sible to have a reliable estimate of their age. We do not find any 
correlation between the disk masses and the stellar ages. The 
first reason for this lack of correlation is that the star evolution 
strongly depends on the stellar mass. Massive stars evolves more 
rapidly than the lower mass ones. We need to separate the stars in 
different ranges of stellar masses to have a meaningful plot. To 
investigate this effect, we have represented different ranges of 
stellar masses with different colors in Fig.|6j3. But even separat- 
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Table 6. Interferometric observations in Herbig Ae/Be stars 



Star 


M^ 


Age 


D 






A 2 


F, 2 


Spec. Index 


M disk 


Inter. 


Ref 




(Mo) 


(Myr) 


(pc) 


(mm) 


(mJy) 


(mm) 


(mJy) 




(Mo) 






ZCMa 


12 




930 


2.7 


3.0(0.7) 


1.4 


18(2) 


2.7(0.9) 


0.12 


PdBI 


1 


MWC 297 


9 




250 


2.6 


79.6(5) 


1.3 


195 (2) 


1.3(0.1) 


0.037 


PdBI, SMA 


1, 2 


RMon 


8 




800 


2.7 


2.1(0.5) 


1.4 


8.7(2) 


2.1(0.9) 


0.009 


PdBI 


1,3,4 


VY Mon 


7 


0.24 


800 


3.4 


1.4(0.5) 


1.4 


46.1(1.4) 


3.9(1.5) 


0.30 


PdBI 


this work 


MWC 614 


5 


0.5 


240 


2.6 


7.6(1.4) 


1.3 


70.8(7) 


3.0(0.8) 


0.04 


OVRO 


5 


BD+61154 


3.9 


1.4 


650 


3.0 


1.7(0.3) 


1.1 


33(1.0) 


3.0(0.6) 


0.17 


PdBI 


this work 


RRTau 


3.6 


1.3 


800 


3.5 


2.1(0.2) 


1.3 


17.9(2) 


2.2(0.4) 


0.19 


PdBI 


this work 


PVCep 


3.5 


0.5 


500 


2.7 


32(1) 


1.3 


280(4) 


3.0(0.1) 


1.14 


CARMA 


6 


VV Ser 


2.7 


3.9 


260 


2.6 


0.7(0.2) 


1.3 


1.4(0.4) 


1.1(0.6) 


0.0012 


PdBI 


7 


AB Aur 


2.6 


3.9 


144 


2.8 


11(2) 


1.4 


85(5) 


2.9(0.7) 


0.03 


PdBI 


8 


HD 34282 


2.5 


5.6 


400 


3.4 


5.0(0.3) 


1.3 


110(10) 


3.2(0.5) 


0.16 


PdBI 


9, 10 


V892 Tau 


2.5 


5.0 


140 


2.7 


39(0.5) 


1.3 


251(10) 


2.5(0.1) 


0.04 


CARMA 


6 


HD 163296 


2.49 


5.79 


122 


2.8 


77.0(2.2) 


1.3 


705(12) 


2.9(0.2) 


0.13 


PdBI 


11 


UXOri 


2.4 


4.2 


450 


2.6 


3.8(0.4) 


1.2 


19.8(2) 


2.1(0.4) 


0.06 


PdBI 


12 


WW Vul 


2.4 


4.2 


550 


2.9 


1.2(0.3) 


1.2 


9.1(1) 


2.3(0.8) 


0.04 


PdBI 


12 


MWC 480 


2.2 


5.98 


130 


2.8 


35.2(0.8) 


1.4 


235(4) 


2.7(0.1) 


0.12 


PdBI 


13, 14 


MWC 758 


2.17 


6.26 


200 


2.6 


6.7(1.3) 


1.3 


56(1) 


3.1(0.6) 


0.04 


PdBI 


15 


CQ Tau 


1.5 


5-10 


100 


3.4 


13.1(0.5) 


1.3 


162(2) 


2.6(0.2) 


0.04 


PdBI 


15, 16 



Refe rences: (l) "Alonso-Albi et al."2009; (2) iManoi et alj|2007l. (3) IFuente et alj|2003l; (4) IFuente et alj 120061; (5) Mannings & Sargent 
(6) iHamidouchei l201d : (7 ) I Alonso-Albi et .all 120081: (8) iPietu et alj|2005l : (9) IPietu et all 120031: (T O) iNatta et all 12004}: (11 ' 
(T2) lNatta et al.ll2001l : (T3) lPietu et alj2006l : (T4) lHamidouche et alj2006l : (T5) IChapillon et alj|2008l: (T6) lBanzatti et alj2011l 



2000 



2007; 



ing the stars in different ranges of stellar masses, we do not find 
any obvious correlation. The first problem is that the number of 
stars in each stellar mass range is very small, specially in the 3— 
7 M range in which we only have 9 objects. Another problem 
is that the estimates of the stellar ages are very uncertain. We 
have discussed in this paper the large uncertainties in the ages 
of BD+61154 and LkHa 198. The uncertainty in the ages of 
some stars is similar to the whole dynamic range of stellar ages. 
Finally, the environment also influences the disk dispersal. Disks 
in clusters seem to be less massive and more scarce than those 
around isolated stars (see e.g. lMann & W illiams 2010). The only 
way to investigate the influence of all these factors in the evo- 
lution of the disks around HAe stars is to have more accurate 
estimates of the stellar ages and observe a large sample of disks 
that allow us to separate the sample in different bins and make 
statistics in each individual bin. 



6.2. Correlation between grain growth and disk masses 

Since it is very difficult to find a correlation between the disk 
properties and the stellar properties, we have tried to investigate 
the disk evolution by correlating two disk parameters, in partic- 
ular the 1.3 mm/2.7 mm spectral index and the disk mass. The 
1.3 mm/2.7 mm spectral index is a grain-gro wth tracer and can 
be understood as an evolutionary indicator. lAcke et ail (2004) 
deduced from the FIR and submillimetre SED that the dust spec- 
tral index de pends on the disk g eometry (flared or flat disks). 
More recently. ISandell et al.l (1201 lb also found a correlation be- 
tween the dust spectral index and the circumstellar mass mea- 
sured on SCUBA maps. Although this result is encouraging, 
it is based on single-dish observations, in which the disk and 
the enveloppe contributions are mixed, and from which the disk 
properties cannot be derived. We have made a compilation of in- 
terferometric observations of disks associated to Herbig Ae/Be 
stars. Only the sources with interferometric measurements at 
both ~ 3 mm and ~ 1.3 mm are considered. In the case of CQ 
Tau, several measurements a re available in the lite rature and we 
used the fluxes measured by Chapillon et al . (2008). More recent 



measurements by Banza ttTet al.l (1201 ll) resolve the disk and are 
not adequate for our statistical study of global properties. 

In Fig. |6j; we show the millimetre spectral index as a func- 
tion of the disk mass (calculated as explained in Sect. 6.1). Disk 
with masses larger than 0. 1 M present a spectral index greater 
than 3, similar to those found in protostellar envelopes. However 
one should keep in mind that in these cases the assumed dust 
opacity of 1 cm 2 g _I at 1.3 mm could be wrong, the real value 
being expected to be closer to ~0.5 cm 2 g _1 . With this value, 
the estimated disk masses would be a factor of 2 larger. Spectral 
indeces lower than 2 are associated with disks with low masses 
(<0.01 M ). In these cases, centimetre-sized grains have grown 
in the midplane and their emissivity is smaller than the canonical 
value. The disk mass could thus be underestimated as well. Disk 
masses in our sample range from 0.003 to 1.14 M , i.e. a factor 
of 380. Taking into account that at most a factor of 10 can be 
due to the unknown grain properties (see Sect 5.5), the range in 
our disk masses would be in anycase larger than 40. Under the 
reasonable assumption that the gradient in the disk mass is due 
to an evolutionary trend with the most evolved objects with the 
less disk mass, this means that only <2.5% of the inital mass in 
dust grains will build up large grains. Another possibility is that 
a significant fraction of the dust mass is locked in large grains 
(> 1 cm) that are not detectable in our observations. 

6.3. Comparison with TTs and Be stars 

There are weal th of evidences of grain growth in disks associ- 
ated with TTs. IRicci et a II (I2010bilal) studied grain growth in a 
large sample of disks in the Taurus and p Ophiucus star form- 
ing regions. They found that the value of f3 is quite constant with 
average value of 0.46. From their results, there is not correlation 
between the dust spectral index and the stellar parameters (mass, 
age). Althougth the stars in the Taurus region are older, the stel- 
lar ages in the p Ophiucus sample are similar to those of the 
stars in our Herbig Ae sample. The characteristics found in the 
disks presented here are also very different from those observed 
in the disks around HBe (> 7 M ) stars (see Paper I). Although 
the stars presented in this paper are older than the HBe stars pre- 
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Fig. 6. a) Disk mass as a function of the stellar mass. The points a re taken from the co mpilation of Paper I, updated with the data 
presented in this work and the PV Cep and V892 data published by Hamidouche (2010). Dashed line corresponds to ~0.04 M , the 
mean value of the disk mass for stars of <3 M . Solid contours indicates stars with masses of 3 M and 7 M . Different colours 
correspond to different ranges of stellar masses: from 1 to 3 M (black), from 3 to 7 M (blue) and > 7 M (red), b) Disk mass as 
a function of the age of the star for the same disks, c) The ~1.3 mm/~3.0 mm spectral index as a function of the disk mass. Only 
disks with interferometric data at ~1.3 mm and ~3.0 mm are considered (Table 6). 



sented in Paper I, their disks have larger masses, and the grains 
in the midplane have not reached centimetre sizes yet. 

Dust coagulation models predi ct that grains gr ow to centime- 
tre sizes in timescales of 10 4 yr dZsom et al.ll2010h . As particles 
grow, the velocity of the particles increases and the fragmen- 
tation becomes more efficient. At some point, the steady state is 
reached and grain growth cease. The maximum size of the grains 
in the steady state depends on the gas and dust physical condi- 
tions (mainly the gas sur face density), the turbulent v elocity and 
the disk temperature (see lBirnstiel et al.1(l2010ll201 ll) for a com- 
prehensive study). The maximum grain size, a max increases with 
the gas surface density and turbulent velocity, but decreases with 
the disk temperature. Disks around HAe stars are warmer than 
those around TTs. This could explain the larger grain sizes in 
TTs. Unless we assume that the turbulent velocity is larger in 
the disks around Be stars, the larger grain sizes in their circum- 
stellar disks is more difficult to explain. 

Another possibility is that the steady state is not adequate for 
the studied stars. Th e stars in our sampl e are very young with a 
high accretion rate. iManoi et al.l (120061) showed that accretion 
in HAe stars decreases substantially at about 3 Myr. Accretion 
could change the grain size distribution resulting in an enrich- 
ment in small grains. Photo-evaporation could also affect the 
grain size distribution, specially in Be stars where the UV ra- 
diation is more intense. Gas and small grains are running away 
from the disk surface at radii larger than the gravitational radius. 
Radial drifts of large particles are not considered in these steady 
state models either. 

7. Summary and conclusions 

We completed our survey of disks around HAeBe stars by imag- 
ing the continuum emission at ~3 mm and ~1.3 mm toward 
BD+61 154, RR Tau, VY Mon and LkHa 198 using the Plateau 
de Bure Interferometer. We have detected the disks associated 
with BD+61 154, RR Tau and VY Mon with disk masses of 
0.035 M Q , 0.05 M and 0.40 M respectively. The disk around 
LkHa 198 has not been detected with an upper limit on the disk 
mass of 0.004 M . The modelling of the SEDs show that the 
grains in the mid-plane of the disks around BD+61 154, RR Tau 



and VY Mon have sizes of ~100, 1000 and 1 //m respectively. 
These massive disks are in the phase prior to the formation of 
protoplanets. Their grains are smaller than those found around 
TTs and young HBe stars. 
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